Octet Spin Fractions and the Proton Spin Problem 
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The relatively small fraction of the spin of the proton carried by its quarks presents a major challenge to 
our understanding of the strong interaction. Traditional efforts to explore this problem have involved new 
and imaginative experiments and QCD based studies of the nucleon. We propose a new approach to the 
problem which exploits recent advances in lattice QCD. In particular, we extract values for the spin carried 
by the quarks in other members of the baryon octet in order to see whether the suppression observed for the 
proton is a general property or depends significantly on the baryon structure. We compare these results with 
the values for the spin fractions calculated within a model that includes the effects of confinement, relativity, 
gluon exchange currents and the meson cloud required by chiral symmetry, finding a very satisfactory level 
of agreement given the precision currently attainable. 



There have been decades of careful experimental investi- 
gation since the original discovery of the so-called proton 
spin crisis by the European Muon Collaboration (EMC) Q]- 
m. The fraction of the spin of the proton carried by its 
quarks currently stands at |5 | 33 =b 3 ± 5 %, if one relies 
on SU(3) symmetry for the octet axial charge, g\. This is 
a dramatic suppression with respect to the value of 100 % 
expected in a naive quark model or even the 65 % expected 
in a relativistic quark model. It increases only marginally, 
to 36 ± 3 ± 5 %, if is reduced by 20 %, as suggested by 
model calculations [6| and a recent lattice simulation [7J. 

A number of possible theoretical explanations have been 
offered, ranging from a key role for the axial anomaly 
Q4J to the effect of gluon exchange currents [[TSUTtII . the 
effects of chiral symmetry [18, 19 ] and, in the light of in- 
sights gained from lattice QCD studies a combination of 
both of these effects {2^. The relatively small values of 
AG, the gluon spin in the proton, found in both fixed tar- 
get and collider experiments lEllEl, have eliminated the 
possibility that the axial anomaly alone might explain the 
observed suppression, although its effect may still be quan- 
titatively significant. 

It is clearly of great interest to find new ways to shed 
light on the origin of this remarkable phenomenon. Con- 
siderable attention is being directed at the measurement 
and interpretation of the generalised parton distributions 
(GPDs) |23|. The moments of these GPDs can be related 
to the quark angular momentum [24J, although there is 
a lively debate over the physical interpretation of those 
moments ll25l[26l . Studies of transverse momentum dis- 
tributions (TMDs) may also offer insight into the orbital 
angular momentum carried by quarks in the nucleon Evll . 
In parallel with these experimental plans, lattice QCD has 
reached a level of sophistication such that one can relatively 
accurately determine the low moments of the GPDs which 
yield the total angular momentum carried by various flavors 
of quarks in the proton |[28j|29l, J^'^'^. The comparison 
of these results with the predictions of quark models after 
QCD evolution (301 appears promising. 



In this Letter we exlore a fascinating new line of inves- 
tigation into the spin puzzle. In particular, we extract the 
quark spin content of the octet baryons from recent lattice 
QCD calculations. These are compared with the predic- 
tions of a relativistic quark model L32.-.34 I which includes 
gluon exchange currents |[T5l [35l [36l and the meson cloud 
required by chiral symmetry |[T8l[33. The variation of the 
suppression of the fraction of the spin carried by quarks 
across the octet is striking and within the relatively large 
uncertainties in the lattice results (which should improve 
significantly in the near future) this is reproduced by the 
model. 

The lattice QCD calculations of the moments of the spin- 
dependent PDFs used here were based upon simulations 
involving 2 + 1 flavors of dynamical quarks, using the 
Symanzik improved gluon action and non-perturbatively 
0{a) improved Wilson fermions 1381 |39l]. These simula- 
tions were performed on a 24^ x 48 volume with lattice 
spacing a = 0.083 =h 3 fm; more details of the specific 
simulation may be found in Refs. ll39ll4Qll . 

Because these simulations were performed at tt and K 
masses between (334,460) and (401,463) MeV respectively, 
it is necessary to extrapolate them to the physical masses. 
This problem of extrapolation has been studied in great 
detail in the literature 141144411 , with the generalization to 
include charge symmetry breaking reported most recently in 
Ref. (451. Our analysis follows that of Ref. (451, where the 
spin dependent moment relevant to the proton spin problem 
is: 

Aqs ^ = f dx{Aq^{x) + Aq^ix)) , (1) 

Jo 

corresponding to the matrix element in the proton of the 
twist-2 operator: 

Ol.^qi.l'^q. (2) 

To lowest order in an SU(3) expansion the matrix element 
of the operator in Eq. ([2]) in a member of the baryon octet 
can be expressed in terms of three coefficients, Aa, A/3 
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and A(t: 



{BmO'^,\B{p)) ^ ^ 
[l\a{BS^B\) + APiBS'^XqB) 
+Aa(BS^B)Tr{Xg)] . 
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Here we have defined 
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with ^ the usual 3x3 matrix constructed from the vr, 77 and 
K pseudo-Goldstone bosons: 
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In addition, the octet baryon tensor, Babc^ is defined through 
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Finally, is Jhe dimensionless spin operator, with 
B-f^YB = -2BS^B, 

The chiral corrections involving vr, 77 and K loops, which 
are non-analytic in the quark mass, are illustrated in Fig.[T] 
These appear at the next order in the formal expansion in 
quark masses. We include both octet and decuplet inter- 
mediate states (the latter indicated by double lines), while 
allowing for the mass difference between them in the nu- 
merical work. As explained in detail in Ref. [45 1, the inclu- 
sion of these loops adds six new 0{mq) fitting parameters, 

An^j*^^ in addition to Aa, A/3 and Act. These are fit 



with the 



to the 24 available lattice data points 
quality of the fit illustrated in Fig.[2j 

Unfortunately, as there are no lattice calculations for the 
A hyperon, we are unable to present results in that case. 
However, for the other members of the octet one can sum the 
values for /S.Uh^ and /\sh to obtain the spin fractions 
carried by the quarks in each octet baryon. Note that because 
the analysis of the renormalization of the lattice operators 
is not yet complete, the absolute values of the spin fractions 
are not known. However, one can compute the ratios of the 
spin fractions for the S and S to that of the nucleon and 
these values are shown in the final column of Table H 

In spite of the fact that at this stage the uncertainties are 
substantial, there is a remarkable degree of variation with 




Figure 1 . Chiral loops included in the present calculation. The 
dashed lines represent the octet of pseudo-Goldstone bosons, the 
lines (double lines) octet (decuplet) baryons and the cross the 
external spin current. 
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Figure 2. Illustration of the fit to the ratio of the lattice moment of 
the doubly represented valence quarks {u in the proton and S^, s 
in the S) and in the E and S hyperons to that in the nucleon. The 
vertical dashed line indicates the physical pion mass. 



the structure of the baryon, with the ratio of spin fractions 
equal to 0.92 ± 0.13 for S : TV, while it is equal to 1.61 ± 
0.33 for S : N. These results clearly do not support the 
hypothesis that the spin suppression observed for the proton 
might be a universal property. It is therefore of considerable 
interest to investigate the predictions of models in which 
the suppression of the spin carried by quarks is dependent 
on structure. With this in mind, we now apply the cloudy 
bag model (CBM), as developed in Refs. |[T5l[T6l[T8l[2Qll. 
to this problem. 

There are three major ingredients of that calculation: i) 
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MIT Bag 


MIT Bag + OGE 


MIT Bag + M. Cloud 


MIT Bag + OGE + M. Cloud 


Model 


Lattice 


N 


65.4 


53.8 


51.9 


43.8 


1.0 


1.0 


A 


77.1 


67.3 


66.4 


58.9 


1.35 (1.33) 






61.5 


50.8 


50.5 


42.6 


0.97 (0.98) 


0.92 (13) 




80.9 


72.3 


72.0 


65.2 


1.49(1.44) 


1.61 (33) 



Table I. Spin fraction (in %) carried by the valence quarks as the corrections discussed in the text are added. The column "Model" 
summarises the full prediction of the ratio of the spin fraction for each hyperon to that of the nucleon (the value in brackets corresponds 
to = 0.8 fm, rather than the default 1 fm). The final column shows the values obtained from our chiral extrapolation of recent lattice 
QCD data for the ratio of the quark spin in the H and E hyperons to that in the nucleon. 
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Figure 3. Gluon exchange current corrections to the spin carried 
by quarks in the octet baryons. The dominant terms are those with 
anti-quark intermediate states. 



the relativistic suppression of the spin of a confined quark 
because of the orbital angular momentum in the lower com- 
ponent of its Dirac wavef unction; ii) color exchange current 
corrections associated with the hyperfine interaction me- 
diated by one-gluon-exchange, and iii) corrections arising 
from the interchange of spin and orbital angular momentum 
when the cloud of pseudo-Goldstone bosons required by chi- 
ral symmetry is included. We briefly outline the calculation 
of each of these corrections within the CRM B2l434iK 

i) Relativity: Within the MIT bag model we take = 
250 MeV as a representative value required to yield the 
observed hyperon masses |49|. Using this value we find the 
spin suppression arising from the lower Dirac component of 
the Is wave function with a bag radius R— 1 fm (0.8 fm) 
to be 0.77 (0.78), compared with the well-known 0.65 for 
a massless quark. Thus the most naive expectation is that 
there should be less spin suppression for the hyperons than 
for the nucleon. 

ii) Gluon exchange current correction: The one- 
gluon-exchange (OGE) force is an essential component of 
spectroscopic studies in most quark models, including the 
MIT bag model. For example, it provides a very natural 
explanation for the A-A^ and S-A mass differences. As 
originally observed by Hogassen and Myhrer I50i it also 
leads to important corrections to spin-dependent observ- 
ables, such as magnetic moments and axial charges, through 
the processes illustrated in Fig. [3] The corresponding cor- 
rection to the nucleon spin was calculated by Myhrer and 
Thomas [15 ], who showed that it reduced the quark spin 
content from 0.65 to around 0.5. We have repeated that 
calculation for the baryon octet. Recall that the dominant 



terms for a spin-dependent external probe are those involv- 
ing an intermediate q — q pair (i.e., with the intermediate 
quark line in Fig. [3] travelling backwards in time). In the 
case of the hyperons there are four different sets of matrix 
elements, labeled fij with {ij) G (g, s) (with q and s a 
light or strange quark, respectively). The second subscript 
refers to the quark hit by the external operator and the first 
to the quark emitting the exchanged gluon. In terms of these 
matrix elements the corrections to the spins of the octet 
baryons, 5X1^, are: 



6^- 



'^fqq 
Vsq 

4 2 4 

~^fqq ~ ~^fsq ~^ ~^fq 



fs. 



3" ^ 3- 
2 4 
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(9) 



where for i? = 1 fm (0.8 fm) our numerical evaluation gives 
= -0.058(-0.058), = -0.049(-0.040), = 
-0.047(-0.049), = -0.039(-0.042). This yields 
the values for the spin fractions shown in the column of 
Table [I| labelled "MIT Bag + OGE". There is only a small 
variation of the size of the correction across the octet, with 
values varying from 12% in the proton to 8% in the S. 

iii) Chiral corrections: For the nucleon the correction 
arising from the pion cloud was first discussed by Schreiber 
and Thomas 1 18 |, who found a reduction of the spin frac- 
tion carried by quarks by 20-30%. At the time there was a 
serious concern about potential double counting if one were 
to combine the OGE and pion corrections, which was only 
resolved (in favor of no significant double counting) after re- 
cent studies of the A-A^ mass splitting in quenched and full 
lattice QCD EOllSll. Already in the late 80's, Kubodera and 
collaborators ITTll combined the OGE corrections with the 
chiral loops for pions, etas and kaons under the assumption 
that there was no double counting problem. 

We have repeated that calculation for the full octet of 
baryons, working strictly within the CBM and using a typi- 
cal bag radius of 1 fm everywhere (including the CBM form 
factors at the meson-baryon vertices). The effect of the 
meson cloud on the MIT bag is shown in column "MIT Bag 
+ M. Cloud" of Table |lj while the final results, including 
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relativity, OGE and the meson cloud corrections, are shown 
in the column labelled "MIT Bag + OGE + M. Cloud". We 
see that once all effects have been included there is a sub- 
stantial variation in the spin fractions carried by the quarks 
across the octet. The meson cloud correction is consider- 
ably smaller in the S than in the nucleon. That, combined 
with the less relativistic motion of the heavier strange quark, 
results in the spin fraction in the S being quite a bit larger 
than in the nucleon. 

The next-to-last column in Table [J shows the ratios of 
the quark spin fractions in the hyperons compared with that 
in the nucleon. That the dependence on the bag radius is 
minimal is illustrated in the second last column of Table |l| 
where the number in brackets shows the model result with 
the bag radius changed from 1 fm to 0.8 fm (the latter 
unrealistic for hyperons) everywhere in the calculation. At 
the current level of precision there is clearly very good 
agreement between the values calculated within the CBM 
and those extracted from lattice QCD. It will be extremely 
interesting to investigate the hyperon spin fractions in other 
models. In addition, this work illustrates the importance 
of further work to reduce the statistical errors of lattice 
QCD simulations of the moments of hyperon spin dependent 
PDFs and to extend them to quark masses closer to the 
physical region. 
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